Design, theory, and measurement of a polarization insensitive absorber for terahertz 

imaging 
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We present the theory, design, and realization of a polarization-insensitive metamaterial absorber 
for terahertz frequencies. We derive geometrical-independent conditions for effective medium ab- 
sorbers in general, and for resonant metamaterials specifically. Our fabricated design reaches and 
absorptivity of 65% at 1.145 Thz. 



INTRODUCTION 

The advent of terahertz (THz) spectroscopy has ush- 
ered in a new field of research as scientists seek to ex- 
ploit the electromagnetic signatures of materials in the 
THz range of the spectrum. Specifically, the ability to 
perform imaging in the THz would have profound im- 
pact on the areas of security [ll |j, [sj, Ij], biology [5|, |g|, 
and chemistry [71, l8[. However, imaging in the terahertz 
is complicated by the a lack of easily accessible electro- 
magnetic responses from naturally occurring materials. 

Electromagnetic metamaterials (MMs) [Hl, Il2 1 are one 
potential solution to overcoming this "THz Gap" 13 1. 
The operational frequency of a given MM design is 
geometrically scalable to other regimes of the electro- 
magnetic spectrum. Various MM structures have been 
shown to naturally couple to either the electric and / 
or magnetic components of light in frequency ranges 
from radio [1'^, microwave 15 1, mm- Wave [16[, THz [13J . 
MIR [17], NIR ;i8j, to the near optical ^. A great deal 
of effort has been exerted to create low-loss metamaterial 
devices such as negative index (NI) structures [15| and 
electromagnetic cloaks[20, l21|. However, the substan- 
tial loss tangent that occurs at the center frequencies of 
metamaterial resonances can be exploited as well. This 
resonant loss phenomenon serves as the starting point 
for the investigation of MMs as narrow-band absorbing 
elements for thermal imaging devices. 

We present the theory, design, fabrication, and mea- 
surement of a single-frequency metamaterial absorbing 
element. The simulated design reaches a peak absorp- 
tivity of 95% at 1.13 THz and the fabricated structure 
reaches a measured absorptivity of 65% at 1.145 THz. 
This absorptivity is comparable to existing MM absorber 
designs [22, l23| , but our design is also polarization insen- 
sitive, which maximizes the absorption of light. 



THEORY 

The absorptivity A(lu) of a given material is given by 
the transmission T{uj) and reflectance R(w) as A{uj) — 
I — T{ll!) — R{uj). In terms of the complex transmis- 



sivity {t) and reflectivity (f), this can be written as 
A = I - |t(w)|2 - |f(w)p. Therefore, A = 1 when 
T = R = 0. In reference [2J], the frequency-dependent 
transmissivity t{Lo) was determined to be dependent on 
the complex index of refraction fi{uj) = ni + in2 and 
impedance Z{lo) = Zi + iZ2 for a slab of length d as: 



i(tj)-i = sin(nfcd) - -{Z + ^) cos{nkd). (1) 

^ Z 

Where k = lo/c and c is the speed of light in vacuum. 
We use the convention where a subscripted 1 and 2 de- 
note the real and imaginary parts of a complex function, 
respectively. 

As Z approaches unity (the free space value), the re- 
flectivity will drop to zero and the transmissivity will be 
determined entirely by n: 



r 



cos{nkd) — isin(nfcd). 



(2) 



Upon substitution of the exponential forms this be- 
comes: 



( — 1 —inikd n2kd 



So the transmission (T — |tp) is 

rp —2n2kd 

Therefore, as ri2 approaches infinity. 



lim T 

712 — *00 



0. 



(3) 



(4) 



(5) 



The combined dielectric and magnetic losses in the sys- 
tem are characterized by 77.2. Therefore the physical in- 
terpretation of the above derivation is that, (in the ab- 
sence of reflections) , the transmission of an electromag- 
netic wave with a given wavevector k through a slab of 
given thickness d is determined entirely by losses in the 
slab. To create a very high absorber it is then necessary 
for Z — I at a point where 712 is large. 



Precise control of n and Z is necessary to realize a high 
absorber. Electromagnetic MMs are prime candidates for 
this task since they can be designed to couple to electric 
and magnetic components of light. This enables precise 
tuning of the complex, frequency-dependent permittivity 
1{io) and permeability /i(w) of a MM slab. The index n 
and impedance Z are in turn given by hiiS) = ■\/e(ti>)/i(a;) 
and Z(ijj) = yj fx{ijj) / i{ijj) . 

MMs are typically highly resonant in e and / or /i, 
where the relevant optical constants approximate the 
form of a complex oscillator in frequency: 



e(w),/i(w) 



, /"o 



F.^^^LO^ 



^0£,M - ^^ - il^ 



(6) 



where F is the oscillator strength, 7 is the damping, loq is 
the center frequency of the oscillator, and £00, /^oo are high 
frequency contributions to e,/i. This form for an oscillator 
describes the frequency response of metamaterials, and 
we term this a "Pendrian" after ref. [12:]. 

We consider a single frequency operation point, defined 
by 0^0, such that when e(a;o) — ft-{^o) then Z[ijJq)=1 and 
"2('-^o) is maximized with a value of: 



"2(^0) 



F' 



(7^0)^ 
which, according to Eq. [H yields 



A{loo) = 1 - Exp(-2:^^i^) 
7 c 



(7) 



(8) 



For the more realistic case, when Eoo 7^ /ioo, then 
Z{ujo) — y/l^oc/^oo 7^ 1 and A is no longer determined 
solely by ^2. However, for large 712, T(wo) remains low 
and R{i-Oo) can be written approximately in terms of ^(0-1) 
[2I such that 



RM = ( 






(9) 



Regardless of Coo and /ioo, an optimal narrow-band ab- 
sorber must maximize F with respect to 7 in both e and 
/x. F is determined by the geometry, filling factor, and 
conductivity of the two metallizations. 7 is determined 
by losses in the metalization and substrate. The optimal 
absorber will then use the MM geometries with the max- 
imum possible filling factor. Furthermore, metallization 
and substrate should be chosen to minimize losses at the 
operation frequency as determined by Eq. [7| above. 



DESIGN 

From a theoretical viewpoint, it would seem simple 
to design a MM that would provide such a response in 




FIG. 1: (color) (a) ERR (b) cross (c) combined ERR and 
cross. The dimensions in microns are a = 84, Li = 52.5, 
L2 = 74, L3 = 64, w= 11, and g= 4. Axes indicate the wave 
polarization and propagation direction. 



e and /i. However, there are several complications to 
the theoretical analysis presented above due to the spe- 
cific properties of MMs. For instance, the periodicity 
inherent to most MMs contributes to spatial dispersion, 
i.e. e = e(w,k) and fi = /x(u;,k).[25[ This spatial dis- 
persion causes the optical parameters to deviate from 
the pure form of Eq. [6l Spatial dispersion also causes 
antiresonances in p, (e) due to resonances in e (/i) [26l |. 
Additionally, conventional electric 21, ^] and magnetic 
[l5| MMs are highly coupled when they share a center 
resonant frequency [29|. Typically, electric metamate- 
rials have higher order electric resonances, and thus in 
general one does not necessarily have the condition that 
£00 = Moo, leading to a nonzero reflectivity as described 
above. In our particular case for the design of a narrow- 
band high absorber, a solution to the latter problem was 
to use two electric resonances; one to raise and lower the 
curve at the operational frequency and one to couple to 
fi to create a region of high absorption. 

The electric responses were provided by a modi- 
fied electrically coupled ring resonator (ERR), shown 
schematically in Fig[TJa). The ERR chosen had four- 
fold rotational symmetry about the pro pag ation axis and 
was therefore polarization insensitive [27|. The lower- 
frequency electric response used to tune the e{Lu) curve 
was driven by the LC loop in the ERR. The higher- 
frequency response used to couple to e(Lu) was created 



by the dipole-like interaction of the metaUizations in ad- 
jacent unit cells. 

The magnetic response was created by combining the 
ERR with a cross structure Fig[ljb) in a parallel struc- 
ture separated by a layer of benzocyclobutane (BCB) 
FiglTJc). The magnetic component H of a TEM wave 
coupled to the center stalks of the two metallizations 
that were perpendicular to the propagation vector, such 
that anti-parallel currents were driven. A polarization- 
sensitive design based on similar parameters has been 
shown at both THz [23| and microwave 22 1 frequencies. 

Practical realization of metamaterials dictates the use 
of a support substrate that complicates the theoretical 
optimization of the structure. This substrate is typically 
thick in the propagation direction compared to the thick- 
ness of the metamaterial and the lateral unit cell dimen- 
sions. This supporting structure introduces ambiguity 
in the definition of the unit cell as well as significant 
asymmetry in the propagation direction, both of which 
complicate the definition of .^e// and ne//.[30| However, 
some portion of the substrate must be included in simu- 
lation because it has a measurable effect on the form of 
the MM resonances due to the dielectric of the substrate, 



31| 



The effect of the substrate on the MM's resonances 
can be shown directly in simulation. A parameter sweep 
of the substrate thickness d shows that the effect of the 
added dielectric on the MM resonance saturates near a 
value oid = ds- For d> dg, the only added effect is a fac- 
tor exp{iy^k{u}){d~ ds) to i(w). As substrates are typi- 
cally chosen for low loss in the frequency range of interest, 
|i(w)| is virtually unchanged and the extended substrate 
only adds phase. The appropriate unit cell boundary in 
the propagation direction is therefore given by dg and 
the MM elements on the substrate. The impedance mis- 
match at the MM-substrate and substrate-air boundaries 
can then be incorporated into the homogenized effective 
medium as these boundaries are within the unit cell. 

The extracted optical constants for the metamaterial 
absorber shown in Fig. [T] are plotted in Fig. [51 The 
real and imaginary components e and fl, as plotted in 
Fig. [21 (a),(c), are complicated by the effects of spatial 
dispersion, as previously mentioned. However, various 
features of the plots enable interpretation of the forms 
of these curves. The lowest frequency feature in €2 at 
Lu = 800 GHz, shown in Fig. [Sjc), is the conventional 
MM ELC resonance. Notice this is accompanied by an 
antiresonance in jlioj), which is characterized by a nega- 
tive imaginary component that peaks at the center fre- 
quency of ei(uj), defined by the cusp near uj — 1.02 THz, 
and in accord with previous work[26j. A second electric 
resonance appears at approximately uj =1.125 THz due 
to the cut-wire response of the cross. This is accompa- 
nied by an antiresonance in jx{uj) centered at the same 
frequency. The resonance in jl{uj) is weak relative to the 
neighboring resonances in e, and thus difficult to observe. 
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FIG. 2: (color) Reflectance (blue), Transmission (red), and 
Absorptivity (green) for the simulated absorber. The vertical 
dashed line indicates the frequency of maximum absorptivity. 
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FIG. 3: (color) Real (blue) and imaginary (red) components 
of e and fx. Real (blue) and imaginary (red) components of 
n and Z. The vertical dashed line indicates the frequency of 
maximum absorptivity. 



However, ^2 has a distinct positive peak at 1.13 THz. 
This, combined with an approximately Pendrian curve 
in /ii, indicates the presence of a magnetic response cen- 
tered a this frequency. There is also a small but distinct 
kink in e{uj) due to the antiresonance caused by the mag- 
netic response. 

Figure [H^b) shows the refractive index n{Lo) and (d) 
impedance Z[uj) and demonstrates how close the struc- 
ture approximates an ideal absorber. At 1.13 THz the 
real impedance is near unity, Zi ^ 1, and the complex 
impedance Z2 is minimized, such that i? ~ 0. As desired, 
the imaginary index, 77.2, is maximized near a value of 



^ 1.75, which minimizes T. This results in a peak ab- 
sorptivity of 95%, plotted as the green solid curve shown 
in Fig. m 



FABRICATION & EXPERIMENT 

The two layer metallization MM absorber sample was 
fabricated on a high resistivity 1 mm thick silicon sub- 
strate. A 3 um thick Si02 layer was deposited on the Si 
substrate using plasma enhanced chemical vapor deposi- 
tion (PECVD). The bottom, cross-shaped metallization 
(30 nm Ti/40 nm Pt/200 nm Au) was patterned using 
standard negative lithography, metal evaporation, and 
metal lift-off. The BCB dielectric (Cyclotene 3022-46, 
Dow) was deposited using two consecutive spin-coat de- 
positions and soft cures in a vacuum oven, resulting in 
a final thickness of approximately six microns. The top 
ELC metallization layer structure and patterning used 
the same process as that of the cross. The two unit cells 
(one for each layer) are shown as the insets to the bottom 
panel of Fig. 4. 

The MM absorber sample was examined experimen- 
tally using a Fourier Transform Infrared (FTIR) spec- 
trometer. Polarized light from a mercury arc-lamp was 
transmitted and reflected from both the sample and a 
reference substrate, and then focused on the detector, a 
liquid helium cooled Si-bolometer. Measurements of the 
sample for both polarizations where characterized and we 
found no deviations within experimental error. Measured 
T{u}) and R{w) were used to calculate the experimental 
A{co). ... 

Fabrication tolerances in the structures resulted in 
some deviation from the theoretically optimized case, and 
are shown in Figure 4. The corners of the ELC and cross 
structures are slightly rounded, and the BCB thickness 
was 6 um rather than the optimal 5.8 um. All of these fac- 
tors were incorporated into the computer model and the 
re-simulated transmission, reflection, and absorbtion for 
the metamaterial absorber are shown in Figure 4, with 
the experimental measurement of T{uj) and i?(w). In 
comparison to the theoretically optimized structure re- 
sults shown in Figure 2, the fabricated sample exhibits a 
shift in the MM resonances. The first transmissive min- 
imum has shifted by lOGHz, while the second minimum 
has shifted by 54GHz, both to higher frequencies. The 
minimum in reflectance has shifted from 1.128 THz to 
1.15 THz. This is expected as electric responses are sen- 
sitive to metallic structure rounding j2l| and the thick- 
ness of the BCB layer partially determines the magnetic 
coupling. However, this change in geometrical parame- 
ters does not fully account for the disagreement between 
Figures 2 and 4. 

At 0.5 THz, T{uj) and i?(w) are well matched in ex- 
periment and simulation, but the forms deviate as the 
curves approach the MM resonances. Specifically, the 



experimental reflectivity reaches a minimum of only 18% 
compared to the simulated value of 2%. Likewise, the 
first and second transmissive minima reach values of 3% 
and 3%, respectively, as opposed to less than 0.1% in 
simulation. Several mechanisms may be behind this dis- 
parity. 

The off-resonant agreement indicates that the discrep- 
ancies are related to resonant forms of the constitutive 
parameters e and p,. As indicated by Eq. [7] lowering F 
or increasing 7 of the resonances translates into a higher 
minimum in T(io) as the peak value of the effective loss 
(^2) is decreased. This also changes the form of Z{u!), 
and therefore the form of R{u!). 

For both the electric and magnetic resonances, 7 is 
primarily determined by the loss in the BCB substrate 
between the ELC and cross. Previous work has deter- 
mined that dielectric losses are the primary mechanism, 
and may be an order of magnitude greater than Ohmic 
loss[22|. Also, because BCB is the primary dielectric that 
tunes the capacitance of each structure, any deviation 
from the nominal value of e =2.5 will lead to variation in 
absorbance. Further, it is well known that the dielectric 
value of many polymer compound have significant depen- 
dence on frequency, especially within the THz range 32 1 . 
As previously mentioned, the "strength of the oscilla- 
tor" F is determined by the geometry, filling fraction, 
and conductivity of the two metallizations. Therefore, 
the deviation of experimental curves in Fig. [3] may be 
caused by a combination of increased loss in the BCB 
and lowered conductivity in the metallizations. 

The best fit of simulation to experimental data is 
shown in Figures 4 and 5. The loss tangent of BCB 
was found to be approximately one order of magnitude 
greater than the nominal value, while the conductivity 
of gold was found to have decreased an order of magni- 
tude. As a consequence of the non-uniform shift of both 
the magnetic and electric resonances, as well as increased 
damping due to loss, n and Z have deviated from their 
optimum values, resulting in a peak absorptivity of 74% 

(Fig. m- 



DISCUSSION 

We now discuss the potential use of the metamate- 
rial absorber as a thermal imager for the THz frequency 
range. Compared to existing THz absorbers, our device 
is narrow band. This enables spectrally-selective appli- 
cations, such as in the detection of ex plos ives. However, 
unlike previous MM absorber designs [23j, our device is 
polarization insensitive. This may be ideal for certain ap- 
plications, as it maximizes absorption for arbitrarily po- 
larized or incoherent light. As a MM device, our design 
is geometrically scalable to different frequency ranges. 
This scalability is limited only by limitations in fabrica- 
tion and loss in constituent materials. The limitations 
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FIG. 5: (color) Experimental (blue) and simulated (red) Ab- 
sorptivity curves. 



design. The theory presented here and the design specif- 
ically show great promise to creating absorbers at any 
decade of frequency. 
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FIG. 4: (color) Reflectance and transmission results for ex- 
periment (blue) and simulation (red) 



of the narrow band, resonant design could be overcome 
by using multiple distinct unit cells [33] or by incorporat- 
ing tunable or frequency agile metamaterial components 

It should be noted that designs presented here are 
bianisotropic - a result of asymmetry in the propagation 
direction - and belong to Shonflies point group C4. 35| 
We have performed simulations, (not shown), in order to 
elucidate the impact of the bianisotropy on absorptive 
properties of the metamaterial. We studied the cross po- 
larization in transmission as a function of frequency for 
the design shown in Fig. [TJ Computer simulations indi- 
cate that the cross polarization is small and achieves a 
maximum of only 10""* over the frequency range of in- 
terest. We also investigated the cross polarizations for 
designs presented in Refs. 2^, l23| and found similar re- 
sults. Thus, for normal incidence radiation, the effect of 
bianisotropy is negligible. 

In conclusion, we have derived general conditions to 
create an absorber based on effective medium theory, and 
for the specific case of MM elements. We have shown that 
such a design can reach absorptivities approaching unity 
within a narrow band. We have also successfully im- 
plemented this approach with a THz frequency absorber 
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